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Finite Element Prediction of 
Structural Failure of Space Payloads 
under Random Excitation using 
Different Solvers 


Abstract— Space payloads are analyzed for their natural frequencies, and for survival under 
inertial loads and random excitation. In this paper, response of a typical space payload under 
Random Vibration is studied using Hypermesh/Optistruct and Nastran/Patran. Quadrilateral 
and Hexagonal elements used in the mesh to obtain higher accuracy and conservative results. 
Gaussian distribution with Kurtosis control 7 is used predict possible failure of the payload. 
Comparison of the results obtained by two solvers is made. A possible explanation to the 
difference in the results is also presented. The explanation helps to approximate results closer to 
theoretical values and the Gaussian distribution of higher order helps to predict failure of space 
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payloads with higher confidence. 


Index Terms— Random Vibration, Dynamic Simulation, Spacecraft payload. 


I. INTRODUCTION 


The space components are required to be designed 
and analyzed with extreme care because they are 
irreparable and need to be maintenance-free. It has to 
withstand with loads caused due to transportation, lift-off 
loads, rocket-motor ignition overpressure, acoustics 
loads, pogo instability vibration loads etc. The Space 
components are usually analyzed and tested for natural 
frequencies to avoid resonance.[ 1] Random response 
analysis is also carried out so that effect of actual 
vibration excitation prevailing during launch can _ be 
simulated. The components are also tested for shock due 
to pyro-shocks and deployment of devises like solar 
panel, antenna etc. Components are subjected to vibration 
through their connection to satellite interface panel, while 
acoustics vibrations are directly transmitted to different 
surfaces of the components. Broadband random 
vibrations are produced due to engine functioning, 
structural response to the broadband acoustic loads and 
aerodynamic turbulent boundary layers. Failure modes 
under consideration in this work have been described in 
Vibration Analysis for Electronics Equipments by 
Steinberg .[2] Possible failure mechanisms displayed by 
components under these loads are described by Stephen 
H. Crandall & William D. Mark in Random Vibrations in 
Mechanical Systems. [3] 

In the present study, a Crystal Oscillator assembly is 
analyzed for structural survivability using normal mode 
analysis, quasi static analysis and random analysis. A 
Crystal Oscillator generates electrical signals at particular 
frequency for Geo-stationary Satellite. Results are 
calculated for the same component with the help of two 
different CAE software Hypermesh/ Optistruct and 
Patran/ Nastran. FEM Model is prepared on_ the 
Hypermesh with the help of 2D shell elements and 3D 
Hex elements. Different types of cards required for 
random response analysis were also studied. 
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Il. RANDOM VIBRATION 


Random Vibration is one of the most realistic kinds of 
excitation in which many frequencies present within a 
bandwidth act simultaneously at any time of excitation 
period. Evaluation of space sub-system under random 
vibration is carried out for simulation of the acoustic 
noise during launch, workmanship levels, pyro-shocks 
etc. The intensity of random vibration loads is given by 
Power Spectral Density (PSD) distribution over the 
excitation bandwidth. [4] 

Methods to Carry out Random Vibration Analysis are: 

a) Direct Integration 

b) Mode Superposition 

In case of direct Integration, the stiffness matrix is 
non-diagonal so there are numbers of equations required 
to be solved. If the structure is large, then it will require 
considerable amount of solver time and computer 
resources. 

In modal superposition, using the orthogonality 
properties of the eigenvector, the stiffness matrix 1s made 
diagonal. Also, the damping matrix is assumed the 
combination of mass matrix and stiffness matrix for 
diagonalization of the matrix. With ortho-normality of the 
matrix, it is possible to reduce the number of equations 
required to be solved and hence, saves time. It can be 
done even on the considerably humble computer 
resources. [5] 

Possible sources of error in Modal Superposition are: 

1. A higher mode frequency is not approximated 
accurately by FEM model even if very fine mesh is used 
for meshing. 

2. Assumption for diagonalization of damping matrix. 

3. Modes extracted are not enough. 


Il]. CASE STUDY: CRYSTAL OSCILLATOR ASSEMBLY 


The function of the crystal oscillator 1s to produce 
electrical signals of a particular frequency. The figure 
shows the OCXO Assembly which also incorporates the 
EPC and OCXO distribution box. 
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Random Vibration Loads - PSD Inputs are as follow: 
TABLE II 
PARALLEL TO MOUNTING PLANE [7] 


Frequency Range PSD Inputs 


20 to 100 Hz + 3dB/octave 


100 to 700 Hz 0.1 g°/Hz 
700 to 2000 Hz -3 dB/octave 


TABLE III 
PERPENDICULAR TO MOUNTING PLANE [7] 


Frequency Range PSD Inputs 
20 to 100 Hz + 3dB/octave 





100 to 700 Hz 0.28 g°/Hz 
700 to 2000 Hz -6 dB/octave 


ie IV. ANALYSIS RESULTS 





A component can fail due to resonance, high stresses, 
Figure 1: CAD Model of Crystal Oscillator acceleration levels or displacement levels. At the internal 
resonance of the system components like relays and 
crystal oscillator due to acceleration, the components will 
not undergo catastrophic destruction but it will impair the 
a erformance of the package or stop the functioning of the 

MATERIAL PROPERTIES USED IN ANALYSIS [6] Pp Pp & Pp & 


; xoungis Poisson Density oe ; ; 
Component | Material Modulus pa (g/cm) There are two kinds of mechanism involved in failure: 
(E) (GPa) (1) the acceleration goes beyond acceptable level for any 
el instant of time. (2) It goes beyond the certain value for 
Base Plate T6 definite period of time. Both of the things can be 


determined by the Time domain analysis of the random 
PCB FR4 25 0.12 1.9 
ee Se 


Major assumptions taken for the analyses are: 

a. Linear Isotropic material properties are assumed 
under all loading conditions. 

b. Masses of connectors and feed through are assumed 
as lumped masses at appropriate locations. 

c. Mass of components over PCB is assumed to be 
uniformly distributed over the surface of PCB. 

d. Fasteners are not modeled in analysis. All interfaces 
fastening joints are simulated with 1-D rigid links. 


The material properties used for analysis are summarized 
in Table 1. 





Also, if the stress level increases beyond the yield 
strength of the material, it may go under structural 
deformation. It is a kind of mechanical failure which 
should be given importance because it leads to 
catastrophic destruction of the system. If two PCBs are 
mounted very close to each other, large relative 
displacements may cause them to strike each other, which 
might result into electronic failure or mechanical failure 
[2]. In this system, to avoid this only one PCB is mounted 
in each box. 

Electrical Signals may go out of suitable range due to 
change in harness, resonance with oscillator, etc. It is also 
important because if this happens the functionality of the 
system is lost [2]. 

To assess any possible failure, analysis results are 


summarized in Table 4 and 5. 
TABLE 2 
NORMAL MODE ANALYSIS RESULTS 


Frequency | Frequency 
Mode (Hz) (Hz) 
No. mee nele Hypermesh Patran 
Optistruct Nastran 


Bending Mode of EPC 213.6 204.1 


Bending Mode of OCXO 
Distribution Box ie 
Bending Mode of EPC PCB 359.8 345.8 


Bending Mode of EPC PCB 386.6 377.1 
Figure 2: FEM Model of Crystal Oscillator Twisting Mode of EPC PCB 434.4 422.9 
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16:01:37 
TABLE V 00-2000, Accelerations, Tx-Basic, .(NON4AYERED) 


RANDOM RESPONSE ANALYSIS RESULTS 
Hypermesh/ Patran/ 
Optistruct Nastran 
Maximum RMS see sa 
Acceleration Tx 
Maximum RMS 
Acceleration Ty 59.12 60.65 
Maximum RMS 


( X- Direction) ( Control PCB) ( Control PCB) 
36.52 MPa 41.3 MPa 
Maximum Von Miss Stress ( Interface ( Interface 


( Y- Direction) between EPC between EPC 
and Base Plate) and Base Plate) 


Maximum Von Miss Stress 32.8 MPa 69.4 MPa 
( Z- Direction) ( Control PCB ) ( Control PCB ) 


Random Reponses for X-directions are shown below: Fas Ee 

















Frequency 


Figure 2: PSD Response for X-direction using Hypermesh/Optistruct 
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Figure No: 6: Von-Mises Stress Plot(X-Direction) using Hypermesh 
/Optistruct 
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Figure 3: PSD Response for X-direction using Nastran 





Figure No: 7: Von-Mises Stress Plot (X-direction) using Nastran 


V. EVALUATION AND INTERPRETATION OF RESULTS 


A. Normal Modes and Inertial Loads 
The Ist modal frequency obtained is 204.1 Hz. For 
this typical payload requirement was to have frequency 


tN above 100 Hz. Hence the margin of safety 1s 104% 
Figure No: 4: RMS Acceleration in X-Direction using 
Hypermesh/Optistruct 





© 2016 RAME IJAEFEA 43 
Research Association of Masters of Engineering www.rame.org.in 


International Journal of Analytical, Experimental and Finite Element Analysis (IJ AEFEA), Issue. 3, Vol. 3, Dec 2016. 


For static inertial load of 20g in X and Y -direction 
and 30g in Z- direction, results are summarized in table 6. 


TABLE VI 
SYSTEM BEHAVIOR UNDER INERTIAL LOADING 


Yield Strength 
(MPa) 


Margin of 
Safety 


Maximum Stress 
(MPa) 


8.67 ( Control PCB ) 413 [1] 


31.42 ( Interface 
between EPC and 
Base Plate ) 


16.18 ( Interface 
between EPC and 
Base Plate ) 


Since margin of safety is positive in all cases, 
structural failure is not predicted due to resonance or 
under inertial loads. 


B. Analysis of Random Vibration Response in Time 
Domain 

It is important to note that different parameters are 
responsible for failure. In case of absence of such 
information, we have to go for RMS values of the 
quantities like acceleration and von-mises stress to check 
survival characteristics of the structure. RMS _ values 
gives some indication of random vibration response, but 
it is not a complete description of the random process. 
Gaussian distribution of excitation is assumed in general, 
which gives a certainty that peak response has a 
probability of 0.27% for crossing 3 times RMS value. But 
a Gaussian distribution assumption for complete 
description of the random vibration in time domain for 
space payloads may lead to inaccurate results. According 
to a study carried out by John van Baren, Gaussian 
distribution with kurtosis control of 7 more accurately 
represents the random process for space payloads. 


Crest factor vs. Kurtosis 
99 994% level 
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Figure No: 8: Crest Factor v/s Kurtosis Control 


In the Gaussian Distribution Process of Kurtosis 
control of 7, Crest Factor is 7.5 which are further used for 
evaluation of results. 


© 2016 RAME IJAEFEA 
Research Association of Masters of Engineering 





e-ISSN: 2394-5141, p-ISSN: 2394-5133, pp 41-45 


Kurtosis = 7 | 
Kurtosis = 3 | 









Probability Density 











Multiple of Sigma 
Figure 9: Comparison between Kurtosis Control 3 and 7 [8] 


Field data goes beyond +3 times RMS value 1.50% of 
the time for space payloads, which is assumed 0.27% in 
case of Gaussian distribution. [8] Gaussian Distribution 
significantly reduces the time spent near the higher 
instantaneous value of acceleration So, It might lead to 
fatigue failure which is caused by the accelerations that 
go beyond +2 times the RMS value. 

The system is subjected to random excitation for 120 
seconds which means that peak acceleration occurs for 
about 7.2 milliseconds. The acceleration goes beyond +3 
times the RMS value for 1.8 seconds which was assumed 
0.324 seconds in Gaussian distribution process. 


RMS Value 
control 7 ) 


The peak value of the acceleration occurs only for the 
instant of time at some nodal point. Hence, it doesn't lead 
to failure because of the inertial effects. RMS values of 
the acceleration and stress are well within the limit, so the 
system is safe under given input conditions. 


TABLE VII 
RESULTS SUMMARY 


Factor Leading to 
Failure 


Peak Value 
(for Kurtosis 





C. Results obtained from Different Solvers 

Different results are obtained from two solvers 
because modal superposition method is used for the 
solution of the dynamic equilibrium equation, and the 
normal modes obtained by two solvers are different. 


VI. CONCLUSIONS 


Crystal Oscillator is analyzed using two different 
solvers and the differences between simulation results are 
studied. Also, study of failure that might occur due to 
random vibration is also done in time domain and the 
structure is found safe for the same. Gaussian distribution 
with Kurtosis control 7 is used for the time domain study. 
A possible explanation for the difference in results from 
the two solvers is also presented. The explanation helps 
to approximate results closer to theoretical values and the 
Gaussian distribution of higher order helps to predict 
failure of space payloads with higher confidence. 
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